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Abstract 


In  the  original  setup  for  the  NURDC  Tower  Experiment,  the 
thermistor  probes  used  to  monitor  the  thermocline  had  time  constants 
of  the  order  of  two  minutes.  Resolution  of  the  data  to  within  0.  5  minute 
was  desired.  Several  wavjr  of  poBt  facto  equalizing  the  thermistor  data 
are  presented  with  the  attendant  advantages  and  disadvantages  of  each. 


1.0  Introduction 


In  the  NURDC  Tower  Experiment,  thermistor  probes  were  used  to 
monitor  the  behavior  of  the  thermocline.  The  time  constants  of  the  probes 
were  to  be  about  0.  5  minutes,  but  an  extra  layer  of  insulation  resulted  in 
time  constants  of  the  order  of  2  minutes. 

There  are  several  ways  in  which  the  response  of  the  thermistor 
probe  can  be  post  facto  equalized  so  that  the  data  will  appear  as  if  the  correct 
thermistor  response  were  present.  The  penalty  for  post  facto  equalization 
will  be  negligible  of  only  phase  equalization  (time  delay)  is  used,  whereas 
there  will  be  an  increase  (about  5  dB)  in  noise  level  if  amplitude  equalization 
is  used  as  well. 

The  crudest  form  of  phase  equalization  is  simply  a  time  shift.  This 
may  be  sufficient.  A  more  sophisticated  approach  is  to  equalize  the  phase 
response  over  the  frequencies  of  interest.  This  can  be  done  in  either  the 
time  domain  as  a  running  weighted  sum  or  in  the  frequency  domain  using  the 
Fast  Fourier  Transform.  The  latter  is  preferred  in  that  the  impulsive 
response  need  not  be  calculated. 

Another  advantage  of  using  the  transform  domain  is  that  the  amplitude 
can  be  equalized  as  well.  This  means  that  not  only  can  the  transients  be 
properly  aligned  but  that  they  can  be  sharpened  as  well  so  that  the  events  can 
be  more  finely  defined  as  to  the  time  of  occurrence. 
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2.  0  Determination  of  the  Thermistor  Response 


In  order  to  equalize  the  response  of  the  thermistor  it  is  necessary 
to  determine  the  response  either  by  a  sec  of  measurements  or  by  a  theoretical 
approach.  The  actual  response  will,  of  course,  be  quite  complicated  and 
reasonable  approximations  must  be  used  to  obtain  a  feasible  model. 


2.  1  Theoretical  Response 

As  a  basis  for  modeling  the  probe,  the  theoretical  response  will  be 
obtained  by  considering  the  thermistor  as  a  sphere  of  insulation  material 
with  the  semiconductor  being  a  small  dot  at  the  center.  The  temperature 
distribution  within  the  probe,  neglecting  specific  heat  effects,  is  then 
governed  by  the  classical  heat  equation  for  spherical  symmetry 


«.  a*  frV-O 


in  which  r  is  the  radial  distance 
/A  is  the  temperature. 


The  initial  conditions  are 

=  ’P(r')  €  t  -  O 

and  the  boundary  conditions  are 

-  7T  e  rm  * 

in  which  ^  is  the  radius  of  the  sphere. 
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The  general  solution 
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When  ^  is  small,  the  infinite  series  becomes 

(ilL  ^  >4rr  =  2Lr 
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and 


^  s  T]  -f  g#  ^77 "7 ]*)  as  expected. 
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When  the  thermistor  bead  Is  small  such  that  ^'^2 ,  the  first 
exponential  term  quickly  dominates  the  remaining  terms  of  the  series  and 

-  £i5 L2* 

4?* 

~r,i- 

Tr  & 


which  for  r"-*  O  becomes 


^  *  ~r,+  3(z-t;)  <2 


-a»*n 


showing  that  the  response,  except  for  very  small  ^  ,  is  governed  by  a 
single  time  constant. 

To  obtain  a  simple  single  time  constant  fit,  the  following  modification 
will  be  made 

/?  ^ 

^  ■  t,+  (7;-t)6 
*  7 7+  Cz-T)  £ 


Thus  the  response  can  be  approximated  by 

-  *l2r\  / 
* 

^  =  ~r,  +  (z-t)  <2 


in  which  the  new  time  axis,  "6*  ,  differs  from  the  old  time  axis,  ^  ,  by  an 


amount 
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The  corresponding  thermistor  response  in  the  transform  domain 


is 

f  _  . 

C  l  +  ST,) 

plus  a  delay  of  2.  .  Formally,  the  transfer  function  is 

-STeJ^a, 

J? _ 

(l  +  ST0) 

An  idea  of  the  accuracy  of  the  approximation  can  be  obtained  by 
considering  the  50  percent  point  in  the  step  transient  response.  The  next 
term  in  the  series  is  then  1/2  (0.  5)"^  or  0.  C3,  The  error  is  thus  about 
6  percent. 

It  remains  to  adjust  the  mcdsl  somewhat  to  account  for  the  actual 
ratio  of  thermistor  vo  insulation  material.  This  corresponds  to  measuring 
temperature  not  at  T  =0  but  at  some  non  zero  radial  location.  The  effect 
is  to  modify  the  delay  by 

~ro  JLi.  Vf 

r< 

Since  this  factor  is  not  known  and  in  actual  £ac .  the  at  delay  could  be 
negligible,  the  response  will  be  taken  as 


J  +  SZ, 


with  the  delay^TJ  handled  separately  should  it  not  be  negligible. 
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2.  2  Experimental  Determination 


The  step  response  of  the  thermistor  probe  can  be  measured  by 
maintaining  the  probe  in  .ir  at  a  temperature  %  until  equilibrium  is 
reached  and  then  quickly  plunging  the  probe  into  water  at  a  temperature  1 r, 


can  be  fitted  to  the  curve  of  temperature  vs  time  by  plotting  (  "7/  ) 

vs  't  on  semi  log  paper  and  fitting  a  straight  line.  The  delay  T(  may 
prove  to  be  negligible. 
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3*  0  Data  Processing  Model  for  the  Thermistor 


Several  models,  of  varying  sophistication,  can  be  used  for  equalizing 
the  thermistor.  The  simplest  is  a  delay,  next  a  phase  equalization,  and 
lastly  both  an  amplitude  and  phase  equalization. 

3.  1  Delay  Model 

The  thermistor  response  can  be  simply  approximated  by  a  single 

delay.  Mathematically,  for  S  (frequency)  small  the  transfer  function  is 
being  approximated  by 

-sr, 

e. 

()  +  ST.) 

This  is  quite  valid  if  most  of  the  energy  or  power  content  of  the  signal  is 
below  the  3  dB  point  of  the  response  since  the  phase  error  at  this  frequency 
is  only  TT* /4  or  i.  e.  ,  the  cosine  of  the  error  angle  is  0.  707. 

3.  2  Phase  Equalization 


-ST. 


eSV» 


e 
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Equalization  of  the  phase  response  can  be  accomplished  in  either 
the  time  domain  or  the  frequency  domain.  The  frequency  domain  is 
conceptually  easier  because  the  impulsive  responses  of  the  equalization 


circuits  do  not  have  to  ba  evaluated. 

The  correct  phase  equalization  is  a  delay  Z} 

angle  Z 

G(oj')  -  u+r~'u>r0 


if  needed,  plus  the 
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This  can  be  done  by  using  the  FFT  to  obtain  the  Fourier  Transform, 
applying  the  delay  and  angle  equalization,  and  then  using  the  Inverse  FFT. 
There  will  be  a  truncation  effect  which  will  destroy  the  usefulness  of  the 
sample  points  near  the  ends  of  each  batch  of  data.  This  truncation  effect 
will  propagate  into  the  data  samples  about  a  time  constant  or  a  minute. 

To  avoid  this,  the  end  points  should  be  discarded  and  the  FFT  time-inter¬ 
leaved  to  cover  this  gap. 

3.  3  Amplitude  and  Pb^se  Equalization 

There  is  very  little  additional  effort  required  to  equalize  both 
amplitude  and  phase  vs  phase  alone.  The  penalty  for  post  facto  equalization 
of  amplitude  is  an  increase  in  the  noise  level.  Indeed,  this  is  what  prevents 
complete  equalization  of  the  amplitude  to  a  flat  response. 

At  this  point  there  are  options.  One  can  equalize  the  2  minute 
response  to  look  like  a  0.  5  minute  response  or  one  can,  within  limits  . 
determined  by  the  deterioration  of  the  noise  level,  equalize  to  other 
responses  or  even  other  bandshapes  than  the  single  time  constant  of  a 
physical  thermistor. 

These  options  bring  up  the  sub,  »nt  of  what  responses  have  good 
transients.  For  example,  if  one  were  to  equalize  to  a  rectangular  low- pass 
transfer  ^unction  with  linear  phase,  the  Gibbs  phenomena  would  show  a 
good  rise  time  but  an  8  percent  overshoot.  Other  bandshapes  would  have 
poorer  rise  t'Tie  but  less  overshoot.  The  single-time  constant  response 
of  a  physical  thermistor  would  have  a  fair  rise  time  with  no  overshoot. 
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In  the  interests  of  credibility  in  the  resulting  data,  only  the  case 

of  equalizing  to  a  0.  5  minute  response  will  be  considered  here.  The 
procedure  is 

1.  Sample  as  long  a  duration  as  the  computer 
can  handle. 

2.  Take  the  FFT. 

3.  Apply  the  transfer  function 

i-h^^r0 
I  +  ^  wZj 

to  each  frequency  of  the  transformed  data 

T„  =  two  minutes  =  120  seconds 
Tj-  half  minute  =  30  seconds 

4.  Take  the  Inverse  FFT 

5.  Discard  the  two  minutes  of  data  a»  each  end 
of  the  data  block. 

6.  Repeat  with  the  new  data  block  overlapping 
the  gaps. 

7.  Add  .a  time  shic't,  if  required. 

The  sampling  rate  should  be  based  on  the  wider  bandwidth  and  signal 
content  to  avoid  aliasing.  As  an  example,  suppose  about  one  percent  aliasing 
is  allowed  based  only  on  bandwidth.  The  corresponding  frequency  is  given 

b>r  wr.'a  /  0  o 

t  *  -i°2.  *  o.frM, 

£rr(3o)  ^ 
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The  sampling  rate  is  twice  this  or  about  1.  0  Hz  or  every  second. 

With  a  4096  point  FFT,  about  4096  seconds  of  data  «*re  processed 
each  time,  i.  e.  ,  68  minutes.  After  the  transfer  functi  tn  equalization  is 
applied  and  the  Inverse  FFT  taken,  the  end  points  must  be  discarded.  The 
extent  of  the  end  points  is  determined  by  the  longer  time  constants  —  about 
2  minutes  at  each  end.  Removing  these  points  still  leaves  an  hour  of  data. 

Processing  time  will  mainly  depend  on  the  I/O  as  well  as  the  time 


to  take  the  FFT  twice. 


The  noi<»e  level  associated  with  this  amplitude  equalization  changes 


by  the  ratio 


^  _ 
1+  (*T/)' 


With 


X  *  u) 


J  '  WTo 

the  ratio  becomes 


1  ^  ( wty)1 


t  +  X  *- 


Ik 

T* 


This  is  a  ratio  of  120/30  or  4  :  1  in  power  -  about  6  dB.  Allowing  that  some 
extra  signal  power  is  present  in  the  wider  bandwidth  (  i  dB),  the  noise 
level  will  probably  deteriorate  about  5  dB. 
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